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A simple approach to the many-polaron problem for both weak and inter- 
mediate electron-phonon coupling and valid for densities much smaller than 
those typical of metals is presented. Within the model the total energy, the 
collective modes and the single-particle properties are studied and compared 
with the available theories. It is shown the occurrence of a charge density 
wave instability in the intermediate coupling regime. 
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I. INTRODUCTION 



The formation of large polarons and bipolarons in polar materials due to Frohlich in- 
teraction with longitudinal optical (LO) phonons has been studied quite extensively since 
the pioneering work of Landau!, Pekari and Frohlichll and now it can be considered a well 
understood problem!. However, a large amount of work has been devoted to the simpler sin- 
gle polaron and bipolaron problems neglecting the effect of the polaron-polaron interaction. 
These effects, on the other hand, are expected to play an important role in heavily doped 
semiconductors! and in many doped perovskites including high T c superconductors!'! that 
are both characterized by large polar effects and densities much smaller than those typical 
of metals. In particular, in the case of superconductors infrared absorption measurements 
suggest the existence of polarons whose properties depend strongly on doping!. 

Recently, the regime characterized by low charge carrier density and strong e-ph inter- 
action has been analyzed. In this regime the formation of a polaron Wigner crystal state is 
favoured with respect to the metallic phase. The stability of such a state in an ionic dielec- 
tric has been studied by the path integral technique pointing out the competition between 
the dissociation of the polarons at the insulator-metal transition and the melting towards 
a polaron liquid!. Moreover two different approaches are known in literature for the metal- 
lic phase of a large polaron system. A perturbative approach, mainly due to MahaniS, 
makes use of the random phase approximation (R.P.A.) in order to obtain a retarded effec- 
tive electron-electron (e-e) interaction due to the exchange of LO phonons0. A variational 
approach, due to DevreeseS, is based on an extension of the Lee, Low and Pines (LLP)§ 
canonical transformation to the many body problem and is able to give the total energy of 
the system in terms of the electronic static structure factor. However, both approaches are 
satisfactory only when electron-phonon (e-ph) effects are weak. 

In order to study the many polaron effects in the weak and intermediate e-ph interaction 
(a < 7) and for polaron densities typical of the metallic phase of the doped semiconductors 
and perovskites we introduce a simple model in which the electron-phonon and the electron- 
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electron interactions are not taken into account at the same level. First the hamiltonian of 
two electrons interacting with each other through Coulomb repulsion and with longitudinal 
optical phonons is introduced and an effective self-consistent electron-electron potential, 
due to LO phonons exchange, is determined within a variational approach. Then the many 
body effects are taken into account considering many electrons interacting with each other 
through the obtained effective potential. Within the proposed model we show that, for 
weak e-ph interaction, many of the known results can be recovered, while, for intermediate 
e-ph interaction, we present evidences for a charge density wave instability as suggested by 
several authorsB'S 

The paper is organized in the following way. In Section II the model is introduced 
and its validity is discussed. The perturbative and the variational approaches are also 
briefly reviewed. In Section III the total energy, the collective excitations and the quasi- 
particle self-energy are calculated for weak e-ph coupling and are compared with the results 
obtained within the variational and perturbative approaches. In Sec. IV the region of 
intermediate e-ph interaction is analyzed and numerical results for the single-particle spectral 
weight function, the collective excitation spectrum and the renormalization coefficient are 
presented just above the appearance of the instability signaled by the complete softening of 
the collective mode. 



II. THEORETICAL FRAMEWORK 

In systems characterized by charge carriers coupled with a polar lattice one has to con- 
sider the electron-electron interaction and the electron-phonon interaction. The Hamilto- 
nian, which can be used to describe these materials, has the following formi: 

H = H + H e _ e + H e ^ ph (I) 

where 
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In Eq. (]l|) the first term represents the kinetic energy of the electrons with effective band mass 
m and the energy of the free longitudinal optical phonons; the second term is the Hamiltonian 
of the electrons interacting with the Coulomb potential screened by the background high 
frequency dielectric constant eoo and the last term describes the electron-phonon interaction 
whose strength, in the Frohlich scheme^, is given by: 

R l/2 /4tt«\ 1 / 2 
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where a is the Frohlich coupling constant, R p = (h/ (2mui)) is the polaron radius and 
V is the volume of the system. In H e _ e and in H e _ p h the q = term is omitted in the 
sums: i.e., we suppose that the self-energy of an uniform positive charge is subtracted from 
Eq.(p. The main properties of the metallic phase of a system described by Eq.(P have 
been discussed within two approaches, proposed by Devreese0 and Mahanll0 respectively. 
In order to make clearer the comparison with the model which we will discuss, we start with 
a brief review of the two approaches. 

The former is a variational calculation of the ground state energy which makes use of 
a generalization of the LLP transformational for many polaron systems. The energy per 
particle of the polaron gas is written as: 

E T = E c + E p (3) 

where E c is the contribution due to the coulomb interaction and E p represents the e-ph 
contribution screened by the electrons, which can be expressed in terms of the electronic 
static structure factor S(q): 



E =-y \MIM (4) 

p Z ^hLO l + h 2 q 2 /(2mS(q)) ' 1 ' 

The latter is a perturbative method of treating the coupled e-ph system and it goes 
beyond the lowest order processes by using the R.P.A. for the e-e potential which contains 
the coulomb repulsion and the e-e interaction mediated by a single phonon exchange: 

V(q,u) = v? + M%D°(q,u) (5) 

where D°(q,u) is the free LO phonon propagator: 

D°(q,oo) = : -. (6) 

LO — Ol>i + 11] LO + LO\ + ir\ 

It is worth noting that the use of the R.P.A. describes the screening of both Coulomb e-e and 
Frohlich e-ph interactions, which are treated on the same footing. In the following we will 
refer to these two models making use of the names variational approach and perturbative 
approach, respectively. 

The model, proposed in this paper, is obtained in the following way. 

First we consider the Hamiltonian of two electrons interacting with the longitudinal 
optical phonons via the Frohlich coupling and repelling each other through the Coulomb 
force: 

P 2 p 2 e 2 

+ E + h.c] (7) 

k 

where p%{r) = {e^ +h.c), R, P, r, p are the position and momentum of the centre of mass 
of the pair and of the relative particle, M and \x are the total and reduced masses respectively. 
Then we obtain an effective potential for the two electrons of eq.(7) eliminating, as we will 
review briefly, the phonon degrees of freedom. Finally the many body effects are studied 
considering a system of electrons interacting with this effective potential. 
The Hamiltonian (7) commutes with 



P^P + ^H^ (8) 
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which is the total momentum of the system. The conservation law of the total momentum 
is taken into account through the unitary transformation: 

U = exp [i(Q - £ fcatag) • R] (9) 

k 

where HQ is the eigenvalue of P t . Following LLP0 we choose the variational trial ground 
state@~El for the transformed Hamiltonian Hi = U~ l HU: 

\ip >= Ui(f)\0 > <p(r) (10) 
where |0 > is the vacuum of and the operator U\ is given by 

U^r) = exp [£C/£(rK " /|(^4)] " ( n ) 
The envelope function 92(7^) is chosen to be an hydrogenic-like radial wave function: 

r r2^ 2/3+3 1 1/2 

(pOjf) = SJU. r p e~ ir (12) 

^ V ; [r(2/5 + 3)J v ; 

where T(x) is the Gamma function. The phonon distribution functions ft are determined 
in a self-consistent way from a functional variational procedure. In particular, the Eulero- 
Lagrange equations for the functions fp 

fi 2 h 2 h 2 h 2 V\m\ 2 

- h v2f ' + h ^ - h Vfi ■ = PsMi (13) 

can be solved exactly for fixed values of the parameters (3 and 7: the solutions are expressed 
in terms of the regular and irregular confluent hypergeometric functionsEHll. Then the 
parameters (3 and 7 of the pair envelope function are fixed, in a variational way, by imposing 
the total energy, e t , to be at a minimum. Since e t reduces to that of two free polarons in the 
LLP approximation!!! when the average distance between the particles is much larger than 
the polaron radius (7 = (3 — 0), the bipolaron binding energy is calculated subtracting from 
the total energy that of two LLP free polarons. In particular it has been shown that the 



bipolaron state exists only if the e-ph coupling constant a is greater than a critical value 
a c = 6 and when rj = eo/e^ is smaller than a critical value rj c = O.OlEHIl 

It is worth to note that in the intermediate coupling regime it is not possible to neglect 
the /•£ dependence on r. In this case the phonons follow instantaneously the relative motion 
of the two electrons and, then, the wave function (10) contains, on average, the retardation 
effects of the e-e interaction mediated by the longitudinal optical phonons. 

Within this variational approach it is possible to obtain an effective e-e potential due 
to the exchange of virtual phonons. In fact, the minimization procedure gives rise to two 
coupled differential equations for f? and (p which are solved self-consistently. With the choice 
(12) for the envelope function it is possible, as mentioned, to solve exactly the differential 
equation for fg as a function of the variational parameters j3 and 7. In the centre of mass 
frame the total energy minimization with respect to ip(r), ^4pJ^> ; gives a Schrodinger- 
like equation for the two electrons from which it is possible to define an effective potential 
depending by j3 and 7 known the functions /J^: 
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Typical effective potentials v(r) are shown in Fig.l. They contain a short range attractive 
term and a long range repulsive term screened at large distances by the static dielectric 
constant e and, in the opposite limit, by the background high frequency dielectric constant 
eoo. It is worth to note that the proposed approach can be, in principle, improved if one 
chooses better and better estimate for the effective potential. 

For the polar materials that are characterized by a value of a smaller than a c , there is 
not bipolaron formation and the variational approach gives 7 = (3 — 0. Consequently v(r) 
coincides with the asymptotic effective e-e potential, i.e., with the effective potential of two 
electrons whose relative average distance is much larger than the polaron radius. On the 
other hand, for a > a c the effective potential supports the bipolaron formation. However, 
it is worth to note that bipolaron formation in the limit of vanishing density does not imply 
the existence of bipolarons at finite densities. 



The proposed procedure allows us to eliminate the phonon degrees of freedom from the 
system, simplifying the treatment of many electron effects, and to investigate larger values of 
the electron-phonon coupling constant with respect to the perturbative approach proposed 
by Mahan. This happens since the LLP transformation^ gives rise to phonon corrections to 
the bare e-ph vertex^, corrections which can be neglected according to the Migdal theoremS 
only when the Fermi energy is much larger than hui (normal metals). 

As mentioned before our approach contains two main approximations: a) the e-e and 
e-ph interactions are not treated on the same footing; b) the variational effective potential 
which contains the e-ph effects is not w-dependent as expected, for instance, in the usual 
Eliashberg theory of sup er conduct ivityS. However, these two approximations are not sig- 
nificant limitations for the physical situation which are aimed to describe. For low values 
of the electron density (low doping) and for both weak and intermediate e-ph interactions 
(a < 7) it is clear that the approach is correct since the e-e effects are less important com- 
pared to the e-ph effects and the physics is controlled by the formation of polarons and, if 
it is possible, bipolarons well described by the LLP variational approach that we adopted. 
A mention deserves the strong e-ph coupling regime (a > 7)i. Recently it has been shown!, 
by using the path integral technique, that the Wigner crystallization of polarons is favoured 
with respect to the metallic state up to critical electron density, n c , of the order of n c ~ 10 19 
cm -3 . For n > n c the system of dielectrics polarons cannot form a liquid state at zero 
temperature and the transition to the metallic phase is driven by the polaron dissociation. 
In this work the regime characterised by very strong e-ph interactions will not be analyzed. 

The regime of larger doping but still such that {u^ < ui, uj® being the plasma frequency 
screened by eo) and weak e-ph interaction also does not present any problem; in fact, as we 
will discuss later, our approach is equivalent to the available theories that do not suffer from 
the restrictions due to the approximations a) and b). 

The validity of our approach is less evident for intermediate e-ph interaction and u>p ~ 
In fact, it is well known that the electron corrections to the bare e-ph vertex tend to sup- 
press the effective phonon-mediated e-e interaction when v F q/u ^> 1, q and uj being the 



transferred momentum and energy. As a consequence, our approach should overestimate 
e-ph effects unless uj is very large. However, when the system is close to a charge-density 
waves instability, it has been shown that the vertex corrections due to the e-e repulsion are 
ineffective also in the limit v F q/uj ^> 1, while the phonon vertex corrections are particu- 
larly relevant^. In other words, it is a reasonable approximation to treat first the phonon 
corrections to the e-ph vertex and then to consider the effects due to the presence of many 
electrons. This recovers the validity of our approach. 

Finally we note that at very large values of the carrier density, typical of the ordinary 
metals, the proposed approach is not able to describe the physics correctly since the re- 
tardation effects of the effective phonon-mediated e-e are relevant, electron screening very 
effective and the phonon vertex corrections are negligible according to the Migdal theorem^. 

The model proposed is studied within the R.P.A.H and Hubbard approximation^! at 
T = 0. Within these approximations the effective interaction between the electrons takes 
the form: 

v , \ v(q) = v(q) n ,x 

e//W,w; i-V(q,cj)/[l + f(q)V(q,u>)] e(q,uu) 1 ' 

where V*(q, uj) = v(q)U°(q,uj) and U°(q,uj) is the lowest order proper polarization propaga- 
tor. The function /(g) takes the value in the R.P.A.il and 



m = 5 (16) 

in the Hubbard approximation!!. The R.P.A. excluding exchange diagrams overestimates the 
contribution of the short range terms by around a factor of two. The Hubbard approximation 
takes into account approximatively these terms making an estimate of the contribution of all 
proper polarization insertions with repeated horizontal interaction lines across the fermion 
loop. In the limit of large electron densities the Hubbard approximation coincides with the 
R.P.A. , the dominant contribution to the correlation energy coming from low momentum 
transfers. 
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The calculation of the self-energy to the lowest order in the dynamically screened inter- 
action V e ff(q,uj) gives: 

E(k,u) = — 1-j HdwJdq &$-j, u -uti-p2L (17) 
h{2ir) Jo J e{q,LUi) 

where G°(q, u) is the propagator for non interacting electrons^ 

r (n . A _ Q (g ~ gg) | gKgf ~ g) /l^ 

In the following we will focus our attention on the properties of the normal state: therefore 
the energy shift is chosen so that the chemical potential of the non interacting system 
H = h 2 kp/2m + €q coincides with that obtained by the equation!!!!: 

M=^ + E(*f,A*-*S) (19) 
Eq.(|l~9"D combined with the above expression for \x gives: 

4 = Z(kF, ^f) • (20) 

Finally to evaluate the self-energy we follow the treatment proposed by Quinn and Ferrel0, 
which retains all contributions coming from the continuum of the electron-hole pair states. 



III. WEAK COUPLING 

The aim of this section is to compare a number of properties, obtained within the model 
introduced before, with the results given by the variational and perturbative approaches. 
In particular we will focus our attention on the total energy, the collective excitations and 
the single-particle self-energy: numerical results will be presented for ZnO, which is char- 
acterised by the following parameters^: m = .24m e , = 4, e = 8.15, fkoi = 73.27 meV, 
a = .849. 



A. The total ground state energy 

The total ground state energy per particle takes the following form: 
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(22) 



In Eg. (|2T|) To and E ex are the kinetic and potential energy, respectively, in H.F. ap- 
proximation; E corr is the correlation energy in the Hubbard approximation and r s is the 
radius, in the Bohr units, of a sphere whose volume is equal to the volume per particle. 
In particular in the limit eo = Coo and m = m e eq. fl2~ID provides the total ground state 
energy per particle of a Coulomb gas. In Fig. 2a the contribution of the e-ph interaction to 
the ground state energy in the variational^! and our approaches are plotted as a function 
of the electron density. In the variational approach this contribution is given by E p in Eq. 
(0). Since in our model it is not possible to express the total ground state energy as the 
sum of two independent terms, the e-ph contribution to Et is obtained subtracting from 
Ex the energy per particle of a Coulomb gas screened by the background high frequency 
dielectric constant eoo. In both models, increasing the density, \E P \ decreases, the screening 
of the Coulomb gas on the polaron self-energy becoming bigger and bigger. However, in our 
approach the effect of the electron screening is slightly less effective for low values of the 
doping. Instead, for very large electron densities, our model underestimates the effects of 
the e-ph interaction. In Fig. 2b the same quantity is plotted as a function of a for a fixed 
value of the electron density, n = 10 18 cm -3 . It is evident that the two approaches provide 
the same results when the e-ph coupling constant is very little (a < 1). 
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B. Collective excitation spectrum 



The collective excitation frequencies of the system are determined by the poles of the 
retarded density correlation function, which occur at the solutions fl q — of the equation: 

e R (q,n q - i% ) =0 . (23) 

In the perturbative approach the dielectric function is chosen to be the sum of three 
contributions!: 

e(q, u) = e 00 + 6 ° ~ e °° 2 - v?IJ?(q, u>) (24) 

1 OJ I UJrp 

where ujt is the transverse optical phonon frequency This expression, which is exact within 
the R.P.A., takes into account the screening effects due to the electron gas, the optical 
phonons and the high energy electronic excitations. Eq.([24j) is a quadratic equation for u 2 
in the limit q — > 0: therefore there are always two roots ojf(q) and oj^iq) that are a mixing 
of the phonon and plasmon excitations. There are two opposite limits in which the many- 
polaron effects can be easily treated: they occur when the plasma frequency (the plasma 
frequency screened by e^) is much smaller or much larger than ui. 

In the former case uj\ tends to ui, the electron gas being not able to oscillate as fast as 
the phonons. On the other hand, since the phonons can follow the motion of the electrons 
uj 2 {q = 0) tends to the plasmon frequency screened by Eq. The situation is different when 
is much larger ui. In fact the electrons see the ionic motion as static so that the static 
electron screening can be used. In this case the frequencies U\(q) and u>2(q) have a different 
behaviour depending on the q value. In fact the electron gas has a characteristic screening 
length q~£, q TF being the Thomas- Fermi wave vector, and the electron-electron interaction 
declines rapidly at distances larger than the screening length, thus for q > q TF the screening 
of the electron gas becomes less effective. For q <C q TF the mode u}%(q) tends to io\ while 
at large values of q it goes to oof since the electron gas is not any longer able to screen 
the electron-phonon interaction. The phonons cannot follow the plasma oscillations of the 
electron gas: then the plasma frequency is screened by the dielectric constant e^. 
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In our approach the description of the collective excitation spectrum is less rich. In fact 
since the retardation effects of the LO phonon-mediated e-e interaction have been neglected, 



Tiujp/ for q — > to the roots of an electron gas screened by the background high frequency 
dielectric constant for large values of q. Particularly at low values of the electron density 
this mode coincides with the plasmon like mode in the perturbative approach. We also 
note that, in our approach, it is not possible to divide the dielectric function into different 
contributions due to the e-e interaction and e-ph interaction. 



In Fig. 3 we show the results of the imaginary part of the electron self-energy lmS(k,u) 
in our model, as a function of to in the Hubbard approximation. Im£(A;, u) is zero for all k 
at hui = Ep\ a change of sign must occur at this point since the damping of electrons and 
holes is opposite in sign. The large peak contribution to ImS(A;, u) is due to the excitation 
of the collective mode. Since the real and imaginary parts of £(fc, o>) obey the Kramers- 
Kronigh dispersion relation, there is a corresponding structure in KeE(k,u): it has a finite 
discontinuity at the same point in which Im£(A;, u) shows a peak. The other contribution to 
ImS(/c, uo) arises from the region where e^(k, u) is not zero, €2 (A;, 00) being the imaginary part 
of the dielectric function, and it is due to the creation of particle-hole pairs. It is also possible 
to evaluate the imaginary part of the self-energy in the perturbative approach in which the 
effective potential between the electrons is the sum of a) the screened coulomb interaction 
and b) the screened electron-phonon interaction. The b) contribution is expressed as the 
product of the screened matrix element Mg/(e^(q, u)) and the phonon Green's function 



there is always 




The collective excitation energy goes from 



C. Imaginary part of the self-energy 



D(q,Lu): 



D(q,u) 




(25) 



where 
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e el (q,uj) = l-v™Tl (q,Lu) . (26) 

Mahan has calculated the one-phonon self-energy term from the electron- phonon inter- 
action replacing the phonon Green's function D(q,u) by the unperturbated propagator 
D°(q,uj) and making use of the static dielectric function of Thomas-Fermi. The result is 
that Im£(/c, u) is zero for energies within u>i around the Fermi energy. The onset of Im£(fc, uS) 
for uj > uji corresponds to the fact that the phonon emission processes are not possible until 
the energy is equal or greater than the phonon frequency Moreover at T = 0, if the 
system is in equilibrium, the contribution to the imaginary part of the self-energy from 
the electron-phonon interaction becomes again zero for Tiu < —Tioji — Ep. The onset of 
Im£(fc, uj) at hu = —Tioji — Ep corresponds to the process in which the electron with energy 
—Ep, jumps in the hole having energy — hui — Ep. These results were first obtained by 
Engelsberg and SchrieffeiS who considered a system of electrons coupled with the acousti- 
cal phonons. We have calculated the self-energy within the perturbative approach making 
use of the R.P.A. approximation and the results are presented in Fig. 4: besides the peak, 
due to the excitation of the coupled phonon-plasmon collective mode, the imaginary part of 
the electron self-energy shows a strong variation in the points +hui, —hui and —Tiuji — Ep, 
but the amplitude of the gap coincides with that obtained making use of the Thomas- Fermi 
approximation for the dielectric constant only in the limit ^> uoi. These results show, as 
it is reasonable, that the static approximation provides a good treatment of the electronic 
effects only when the plasma frequency is much larger than the characteristic phonon en- 
ergies. Similar results have been obtained by Das Sarma et alJUl who have calculated the 
quasi-particle damping rate which is proportional to the quantity T(k) = ImS(fc, k 2 ). They 
showed that the static screening is a good approximation at very high electron densities, 
whereas the dynamically screened decay rate lies intermediate between statically screened 
and unscreened results. In Fig. 5 T(k) is plotted as function of k for the ZnO at n = 10 17 
cm -3 (uj^ 3 uji), both in our model and in the perturbative approach making use of the 
R.P.A. . In particular in the region close to the Fermi surface the damping rate goes like 
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(E°(k) -n) where E°(k) is the quasi-particle energy: this is a simple consequence of the 
Pauli principle restrictions. The peaks in T(k) are caused by processes corresponding to the 
excitations of the two renormalized modes u>i(k) and uo2{k). Moreover it is evident that at 
low electron densities uj? 3 <C uji our model neglects, with respect to the perturbative ap- 
proach, excitation processes involving high energy, whereas it provides a good treatment of 
the characteristic polaronic energies. We draw the same conclusions evaluating the polaronic 
corrections to the band edge which, in the Rayleigh-Schro dinger perturbation theory, are 
given by E p = ReSi(A; = 0, u = 0), £i(fc, u) being the one-phonon self-energy term from the 
e-ph interaction. In Fig. 6 we present the results for E p in the perturbative approach making 
use of the R.P.A., static R.P.A. and the Thomas-Fermi approximations. It is clear that the 
Thomas-Fermi approximation overestimates the screening effects at low densities, giving a 
smaller polaronic binding energy than the dynamically screened results. In the opposite 
limit, instead, the two approximations provide the same results. 



D. The spectral weight function 

The spectral weight function is the imaginary part of the electron propagator and there- 
fore it can be expressed in term of the real and imaginary parts of the self energy £ ret (&, c<j)il 

A(k, U) = ^rJM . (27) 

[uj — E®/h — ReT, ret (k, uj)] + [lmS ret (k, u)] 
The Fig. 7a, 7c show A{k,u) in the Hubbard approximation for two different values of i]. 
The energy E(k) and the damping rate T(k) of a quasi-electron (k > kp) and a quasi-hole 
{k < kp) are the poles of the analytic continuation of G(k, to) into the lower right and upper 
left half plane for uj : 

n2 k 2 r,^(, E(k)+iT(k)\ 



and 



m= i^{k,m±my (29) 
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In particular for long lived single particle excitations the energy E(k) and the damping rate 
T(k) are given by: 

S(*) = *£ +BaE (*,M) ,30) 

and 

m = , Im f if ' . (3D 



, __ dRcE(fc,q/> 

1 a 



The numerical solution of Eq. (j30|) shows the existence of three roots. The one at highest 
energy is the regular quasi-particle solution (screened polaron), which is very little shifted 
from the Hartree-Fock value. It describes electrons surrounded by a dynamic polarization 
cloud due to the effective e-e interaction: the self-energy at q = of this quasi-particle 
increases with the particle density. Another solution appears at the finite discontinuity 
of ReS(fc,u;) and it doesn't give contribution to the spectral weight function having the 
imaginary part of the electron self-energy a peak at the same energy. To understand the 
physical meaning of the other solution it is convenient to separate the contribution of the 
exchange term from the self-energy: £(&, u) = H exc (k, u) + £*(&, u). The contribution given 
by S*(fc, a;) can be viewed as that of a particle coupled to the density fluctuations propagator 
S(k,co), whose spectral representation is: 

S(k,u>) = — - f duilm— — tj . (32) 

On the other hand, the single-particle self-energy of a fermion coupled to a boson field 
with the Bose propagator: 

D(k,u)= 2 2 "^) (33) 
uj l — io l [k) + id 



and the effective coupling: 
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takes exactly the form £*(&, ui) if the further approximation of single plasma pole for the 



as the contribution of a particle interacting with a boson field: the plasmon-phonon. The 
quasi-particle formed by an electron surrounded by a cloud of virtual excitations (plasmon- 
phonon) is usually called plasmapolaronil. 

The quantity E — t* 0r) where E is the q = self-energy of the plasmapolaron, i.e. the 
solution of the equation E = Re£(0, E), and e$ T is the shift in the chemical potential for a 
Coulomb gas screened by e^, is plotted as a function of the electron density in Fig. 8. It is in 
good agreement with that calculated starting by the Hamiltonian of an electron interacting 
with the phonon and plasmon boson fields and using a variational approach based on the 
LLP transformation^. The differences between the two plots are due to the presence, in the 
self-energy, of the continuum of the electron-hole pair states, exchange and correlation terms. 
These contributions are not considered in the variational approach proposed in ref. [37]. 

In Fig. 7b, 7d the spectral weight function in the perturbative approach is plotted for 
the same values of i], n and k. It is well known that, in the perturbative approach, A(k,u) 
for a single electron coupled to the longitudinal optical phonons shows, at T = and 
for k — 0, a delta-function peak and an incoherent contribution, with very small spectral 
weight, which starts at uj = cjJI. The effect of the e-e interaction is to split the single 
delta-function peak in two lorentian functions that describe two quasi-particle: the polaron 
and the plasmapolaron. The comparison between Fig. 7a and Fig. 7b points out, again, that 
our approach takes well into account the characteristic polaronic energies and neglects the 
excitation processes involving high energies {uj > Mi). This is a reasonable approximation 
when < uj\. 

It is worth to note that, increasing the charge carrier density, the polaron, an electron 
dressed by a cloud of virtual phonons, evolves towards a new quasiparticle, an electron weakly 
renormalized by the scattering with the phonons, i.e., a quasiparticle with reduced energy 
and mass with respect to the values E p (q = 0) = ~ahio\ and m* = mj (1 — a/6). Instead 
the plasmapolaron evolves towards an electron dressed by a cloud of virtual plasmons. 



density fluctuations 




Therefore 
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IV. INTERMEDIATE COUPLING (a ~ a c ) 



A. Collective excitation spectrum 

It is well known that the retarded dielectric function is analytic in the upper half of the 
complex plane provided!! e R (k, 0) > 0. There are examples of physical systems for which 
the linear response function violates the causal requirement: in this case the linear retarded 
dielectric function can no longer describe the behaviour of the system correctly. It has been 
shown that the temperature dependent correlation function in the R.P.A. for an interacting 
many-particle system, e R (k,u), when the interaction is sufficiently attractive, has forbidden 
zeros on the imaginary axis of the complex frequency planed. These poles of the linear 
response function appear even in the classical limit where they correspond to the transition 
of the system from a gaseous to a liquid state. A similar situation is found in the the case 
of a superconductor!^. While in the perturbative approach the causal requirement is never 
violated, even if the interaction electron-phonon is extremely strong, in our model, if the 
value of the coupling constant a is sufficiently large, the linear response function possesses 
a pair of imaginary poles. It is interesting to note that the same type of instability has 
been suggested by Di Castro et al. for a large class of systems of interacting electrons and 
phonons as due to very ineffective electron corrections to the e-ph vertex0. 

The results which will be shown here and in the following subsections have been obtained 
choosing the model parameters in such a way the system is just above the appearance of 
the instability. The parameters chosen are: m = m e , rj = .01, R p = 10 A, afiuj = .016 Ry, 
(a ~ a c ). 

In Fig. 9 we present the numerical results for the collective excitation spectrum in the 
Hubbard approximation. The collective mode energy softens for a critical wave vector q 
indicating strong correlation between the electrons. If the attractive potential is sufficiently 
strong the collective energy softens completely. This softening, present in a finite range of 
densities 10 19 — 10 21 cm -3 , means that the system becomes unstable with respect to the 
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formation of the charge density waves. Similar results are also obtained in the Hubbard 
model with on site attraction!! and using the slave-boson approach for the infinite U three 
band Hubbard mo delS 



B. The spectral weight function 

The Fig. 10 shows A(k, u>) for k = O.Olkp and for different values of the electron density in 
the Hubbard approximation. Beside the peaks already present in the weak coupling regime 
the Eq.(|30D exhibits a new solution with energy between the plasmapolaron and quasi- 
particle branches. It appears only if the strength of the attractive potential is sufficiently 
large and the damping rate of this excitation increases very quickly as a function of k. 
Then it exists as a well-defined excitation of the system only for moments well inside the 
Fermi surface. As mentioned before this excitation is related to the attractive part of the 
potential and it is present only close to the instability signaled by the complete softening 
of the collective mode energy. This interpretation is also confirmed by the behaviour of the 
renormalization coefficient of the one-electron Green's function: 

Z ( k ) = - ^em ' ( 35 ) 

At k = kp-, Zikp) is the magnitude of the discontinuity at k = kp in the momentum 
distribution n (fc)i. In Fig. 11 Zikp) is plotted as a function of the electron density and it is 
compared with the renormalization coefficient of a Coulomb gas screened by the background 
high frequency dielectric constant In the range of densities in which the new elementary 
excitation appears, Zikp) has a minimum; increasing the electron density it tends to the 
value of a Coulomb gas screened by too. These results indicate that the charge density wave 
instability is a source of a strong quasi-particle scattering and suggest that the new peak in 
the spectral weight function could be a metastable double occupied state (bipolaron). 
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C. The dispersion curves 



We end this section showing the dispersion curves of the three excitations Ei(k) in the 
region of k in which the decay is due only to the excitations of particle-hole pairs (Fig. 12). 
They have been obtained following the poles of the analytic continuation of G(k,u). In 
particular, for k less than /cp and energies less than fi, these poles describe the state obtained 
by creating a hole in the interacting ground state. The curves Ei(k) are drown only for k 
values which do not involve the creation of collective excitations. In fact the decay for the 
quasi-particle band due to the creation of collective excitations starts for a value of the 
electron momentum which satisfies the conservation of the energy: 

*£=*S +M *_*,) (36, 

u(k) being the collective excitation spectrum of the system, while the critical momentum 
which characterizes the decay of a plasmapolaron is given by the condition: 

E(k) - fx < -%uj{k + k F ) . (37) 

Both the quasi-particle and plasmapolaron have approximatively quadratic dispersion laws. 



V. CONCLUSIONS 

A model for an interacting gas of large polarons has been introduced. In the model the 
electrons interact each other with a non retarded effective potential obtained within a vari- 
ational approach!!! valid for large polarons in the intermediate e-ph coupling. The effective 
potential is able to recover the large polaron properties in the limit of zero density. Then the 
many-body problem has been studied within the R.P.A. and Hubbard approximations. The 
ground state energy, the single-particle self-energy and the collective excitation spectrum 
for a system of interacting large polarons have been calculated at densities typical of the 
metallic phase of doped polar semiconductors and doped perovskites. The numerical results 
have been compared with those provided by two traditional approaches in which the e-e 
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and e-ph interactions are taken into account at the same level. It has been shown that the 
polaron self-energy is strongly reduced compared to that obtained in the zero density limit. 
For intermediate values of the e-ph coupling constant the collective mode energy presents a 
softening in a finite range of densities at a finite value of q signaling a CDW instability as 
suggested in recent works about the e-ph interaction in presence of strong correlations^ and 
the single-particle spectrum shows a new elementary excitation, a bipolaronic metastable 
state, with energy between the plasmapolaron and the quasi-particle branches. 
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FIGURE CAPTIONS 

Fig.l. The effective self-consistent electron-electron potential. The energies are given in 
units of hui and are measured from —lahuji. The parameters of Fig. la correspond to 
the characteristic values of ZnO. 

Fig.2. a) Electron-phonon contribution to the ground state energy per particle in the 
variational (dotted line) and in our model (solid line) as a function of the electron 
density. The energies are given in units of Hcuf, b) The same quantity is plotted for a 
fixed value of the electron density, n = 10 18 cm -3 , as a function of the e-ph coupling 
constant. The parameters are the same of ZnO, except eoo, that varies from 1.4 to 
8.15. 

Fig. 3. The imaginary part of the self energy for k = .5kp as a function of the electron 
energy hu measured from at n = 10 17 cm -3 . ImE(^,w) is given in units of 4Ep. 
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Fig.4. The imaginary part of the self energy for k = .5kp as a function of the electron 
energy Tiuj measured from ej at n = 10 17 cm -3 in the perturbative approach. In the 
region of the characteristic polaronic energies these results are very similar to that 
obtained in our model reported in Fig.3. ImE(fc,u;) is given in units of 4Ep. 

Fig. 5. The imaginary part of the quasi-particle self energy Im£(/c, k 2 ) as a function of the 
momentum k of the electron at n — 10 17 cm" 3 in our model (solid line) and in the 
perturbative approach (dotted line). The energies are given in units of fiuii. 

Fig. 6. The contribution of the electron-phonon interaction to the correction of the 
band edge in the perturbative approach making use of the R.P.A.(solid line), static 
R. P. A. (dashed line) and Thomas- Fermi (dotted line) approximations as a function of 
the electron density. The energies are given in units of fiuii. 

Fig. 7. The one-particle spectral function in the range of the characteristic polaronic energies 
for k = .5kp and at n — 10 17 cm -3 , in our model and in the perturbative approach 
for two different values of i] (a, c and b, d respectively). In both cases e$ is of order 
of 6q ~ — ahuji. The energies are given in units of AE F . The parameters of Fig.7a,7c 
correspond to the characteristic values of ZnO. The spectral weights relative to the 
different peaks are indicated. 

Fig. 8. The plasmapolaron self-energy as a function of the electron density in our model 
(solid line) compared with that obtained in ref. [37]. The energies are given in Ry and 
are measured from —ahcui. 

Fig. 9. The curves a) and c) represent the collective excitation spectrum of the Coulomb 
gas screened by e and respectively at different electron densities (n = 10 19 cm -3 , 
n = 10 20 cm" 3 , n = 10 21 cm" 3 ). The particle-hole excitations lie between the dotted 
lines. The curve b) represents the energy of the collective mode as a function of the 
wave vector q using the potential plotted in fig.l. The energies hQ(q) are given in 
units of AEp- 



Fig. 10. The one-particle spectral function for k = .01k F as a function of the electron energy 
at different charge carrier densities (n = 10 19 cm -3 , n = 10 20 cm -3 ). The energies are 
given in units of 4Ep- 

Fig. 11. The renormalization coefficient of the one-electron Green's function at k = kp as 
a function of the density. The solid line represents Z{k F ) for a Coulomb gas screened 
by the background high frequency dielectric constant e^. 

Fig. 12. The dispersion curves at different electron densities (n = 10 20 cm -3 , n = 10 21 cm~ 3 ) 
in the region of k in which the decay is due only to the particle-hole pairs excitations. 
The energies are given in units of 4E F . 
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